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^ Abstract 

(N 

The present density functional theory study addresses the question whether 
the presence of H2O influences the catalytic activity of small gold clusters, 
^ Aui_4/MgO(100), towards the oxidation of carbon monoxide. To this end, we 

4^ studied the (co-) adsorption of H2O and CO/O2 on these gold clusters. The 

^ ground state structures in the presence of all three molecular species, that we 

found, are Aui02/MgO and Au2-4CO/MgO with H2O adsorbed on the surface 
^ in the proximity of the clusters-molecule complex. In this configuration the 

J t catalytic activity of Aui_4/MgO is indifferent to the presence of H2O. We also 

found that a stable, highly activated hydroperoxyl-hydroxyl complex, 02H-- OH, 
^ can be formed on Aui^s/MgO. For the catalytic active system Aug/MgO, it 

has been predicted that this complex opens an alternative catalytic reaction 
pathway towards CO oxidation. Our results suggest that this water mediated 
Q catalytic cycle is unlikely to occur on Aui^a/MgO. In the case of Aui/MgO 

O the cycle is interrupted by the dissociation of the remaining (OH) 2 complex 

after forming the first CO2 molecule. On Aus/MgO the O2H • • OH complex is 
likely to decompose upon the impact of a CO molecule, since three of its bond 
^ dissociation energies are comparable to the reaction barrier of the CO to CO2 

On oxidation. 
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1. Introduction 

Gold, known and treasured for being chemically inert in bulk, received new 
^ attention after the discovery [TJ [2] that supported nanometer-sized gold parti- 

cles have a substantial catalytic activity towards low-temperature oxidation of 
rS carbon monoxide [Hill [5]. Even sub-nanometer gold clusters, consisting of a 

^ few atoms only, showed an interestingly high activity towards CO oxidation, 

both in the gas phase and on metal-oxide supports [3 [3 [H [T^ . During the 
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last decade much, predominantly theoretical, work has been carried out in this 
field, revealing for instance the importance of the cluster-substrate interaction 
[TTl [12] and of the defects in the substrate [I3l [14] . The influence of other chem- 
ical species on the catalytic activity of supported gold clusters has been studied, 
e.g. for chlorine, which hinders the catalytic activity [15 . 

The effect of moisture on the catalytic activity of supported Au nanoclusters 
towards CO oxidation is less conclusive. Both, enhanced activities of up to two 
orders of magnitude for Au/MgO [16 as well as inhibition of the reaction on 
Au/Ti02 [IZ] have been found in experiments. Magnesium oxide is one of the 
most studied substrates used for gold cluster deposition because of its stability 
and simplicity to handle both in theory [T8| and experiment [19] . The catalytic 
activity of small gold clusters, Au^ (n = 1,..., 20), towards CO oxidation has 
been observed [6 . Especially, Aug/MgO gained much attention [9]. So far, 
it is the only system of this type, where the influence of H2O on its catalytic 
activity has been studied. In Ref. [20] it has been predicted that H2O can open 
an alternative reaction pathway for the oxidation of carbon monoxide, hence 
promoting the activity of Aug/MgO. 

Recently, the co-adsorption of CO and O2 on Aui_4/MgO(100) has been sys- 
tematically studied [21 . In that work the absence of catalytic activity towards 
CO oxidation of Aui^2/MgO and the small activity of Au3^4/MgO, observed 
experimentally in [6], could be explained by theory. 

Here we address the question whether the presence of H2O in these model 
systems influences the adsorption behavior of CO and O2 on Aui_4/MgO(100) 
and if a water promoted catalytic activity, as suggested in Ref. [20 , can be found 
even for the smallest gold clusters on a regular MgO terrace. We concentrate our 
affords on Aui_4/MgO, because of their restricted number of isomers, which al- 
lows us to reliably determine the ground state structures. These sub-nanometer 
gold clusters are immobile on a regular MgO terrace at low temperatures due to 
their size-dependent diffusion barriers Ediff of 0.19 to 0.62 eV [22 and adsorp- 
tion energies Eads of -0.89 to -1.75 eV [21]. Our model system is appropriate 
under the assumption of a low defect concentration on the MgO surface. 

The adsorption and diffusion of a wide range of molecular species on regular 
MgO terraces has extensively been studied by theory. The reported adsorption 
energies are -0.08 eV for CO [23 and -0.04 eV [24J for O2, respectively. Since 
CO and O2 bind significantly stronger on Aui_4/MgO than on bare MgO [21 , 
we concentrate on the (co-) adsorption of these molecular species from the gas 
phase. The H2O molecule has been found to bind stronger to MgO and being 
less mobile, i.e. Eads = -0.37 eV and Ediff = 0.13 eV [25l[26]. Its adsorption on 
Aui_4/MgO has not been studied so far. 

The paper is organized as follows. After summarizing the computational 
details in section|2| we will, as a first step, present our results of the adsorption of 
H2O on Aui_4/MgO, section [3] This allows us, in a second step, to investigate 
the possible co-adsorption scenarios of CO or O2 with H2O on Aui_4/MgO, 
section [4| In section |5] we address the question, whether a water promoted CO 
oxidation cycle, as proposed in the case of Au« /MgO,[2Q] is possible even for the 
smallest gold clusters. Finally, section |6] summarizes our results and concludes 
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this work. 



2. Computational details 

The scalar-relativistic ab-initio DFT calculations were performed using the 
projector augmented wave (PAW) [27, 28 method as implemented in VASP [291, 
|3Q] . The exchange-correlation interaction was treated in the generalized gradient 
approximation (GGA) in the parameterization of Perdew, Burke, and Ernzerhof 
(PBE) [31 . A cut-off energy of 600 eV was used and a Gaussian smearing of a = 
0.02 eV for the occupation of the electronic levels. Spin-polarization was taken 
into account for all calculations. The MgO(lOO) surface was modeled in a super- 
cell approach with a two monolayer thick 3x3 MgO slab. The unit cell was 
constructed using the equilibrium lattice parameter of MgO (4.235 A) obtained 
in the corresponding bulk calculations. The repeated slabs were separated from 
each other by 27 A of vacuum. A Monkhorst-Pack F-centered 2 x 2 x 1 k- 
point mesh was used for the structural relaxations. The relaxation cycle was 
stopped when the Hellmann-Feynman forces had become smaller than 5 • 10~^ 
eV/A. The char ges were calculated by means of the Bader analysis [32 . For 
an exothermic adsorption of A on adsorbent S the adsorption energy, £^ads = 
£^o[A/S] - ^o[A] - Eo[S], is negative. 

3. Adsorption of H2O on Aui_4/MgO 

The gold clusters Aui_4 stand on the regular MgO(lOO) terrace and form 
bonds with oxygen atoms Og, which donate 0.3 e~ per Au-Og bond to the 
clusters [El EH [21]. 

Starting from various initial configurations, we identified two possible sets of 
adsorption sites for H2O. First, on the MgO surface in the proximity of Aui_4, 
Fig. [1] and, second, on top of Aui_4/MgO, Fig. [2j The adsorption energies for 
both cases are shown in Fig. [3] (a) together with Eads on bare MgO. 

We found that the binding mechanisms of H2O in the proximity of Aui_4 
and on bare MgO are rather similar: in both cases the water-oxygen atom binds 
to a surface magnesium Mg [34| |35l [25]. It receives a small addition charge 
of up to 0.1 e~ from the surface oxygens surrounding this Mg. The charge 
donations from the surface oxygens to the gold clusters and the water molecule 
are virtually unaffected by their mutual proximity. Also the structural changes 
in the Au clusters are rather small, i.e. bond lengths change by less than 1.5%. 

In the proximity of Aui_4/MgO two effects infiuence the adsorption energy 
of H2O: first, the Coulomb repulsion between the negatively charged cluster and 
molecule and, second, the orientation of the water molecule's electrical dipole 
towards the clusters. At the binding site on top of the next-nearest Mg from the 
gold cluster H2O gains most energy from the combination of these two effects. 

The energy cost of separating H2O from the Au cluster, Egep = Eo[Aun/MgO 
with H2O in proximity] - Eo[Aun/MgO] - Eo[H20/MgO] shown in Fig. |3](b), 
is identical to the difference E^J^^^ - E^^g® of the H2O adsorption energies close 
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and further away from the cluster. Therefore the adsorption energy of the gold 
clusters is not affected by the proximity of the water molecule and the energy 
is gained from the orientation of the electrical dipole, only. 

The other (meta-) stable adsorption sites for H2O are on top of the gold 
clusters, see Fig. [2] and Fig. [sj^a) for the adsorption energies. The adsorption 
on top of the clusters shows pronounced cluster-size effects: depending on the 
parity of the number of Au atoms in the cluster H2O forms an Au-H (odd) or an 
Au-0 bond (even), see Fig. [2] The oscillations of properties due to the changing 
parity of the number of Au atoms were previously observed for low dimensional 
gold structures ^3 • The different binding mechanisms for the top sites are 
readily explained by the electronic structures of Aui_4/MgO: gold clusters with 
an odd number of atoms have an unpaired 6s electron that forms a weak bond 
with hydrogen's Is electron. On the other hand the Sd-states are the highest 
occupied states of the clusters consisting of an even number of atoms, forming 
a stronger bond with the 2p electrons of oxygen. 

These different odd-even bonding mechanisms are, for instance, also reflected 
in the oscillation of the gold-surface oxygen bond length dAu-o^ and the charge 
redistributions to the water molecule. For Aui^a dAu-o^ decreases by 2.5%, 
whereas it increases by 2% for Au2,4. The most remarkable water-induced struc- 
tural change is the breaking of the Au-Au bond between the two top gold atoms 
of Au4 with H2O in the top position. Adsorbed on top of Aui^a up to 0.05 e~ 
are redistributed to the water molecule, while the molecule loses up to 0.06 e~ 
when adsorbed on Au2,4, see the Bader charges in Fig. [2] 

4. Co-adsorption of H2O and CO/O2 on Aui_4/MgO 

Recently the (co-) adsorption of CO and O2 on Aui_4/MgO(100) has been 
systematically studied [21 . Here, based on the analysis of the individual ad- 
sorption energies of H2O, CO, and O2, we can determine the cluster- molecule 
ground state structures in a mixed atmosphere of all three molecular species. 
In Figures [4] and [5] the adsorption energies of CO and O2 on Aui_4/MgO with 
and without of H2O are compiled. 

As mentioned above, CO and O2 adsorb weakly on a regular MgO terrace, 
most favorably on top of Mg, with Eads = -0.09 and -0.03 eV, respectively 
(see dotted horizontal lines in Figs. [4] and [5|. On top of the gold clusters 
Aui_4/MgO, on the other hand, CO and O2 bind stronger (open circles in Figs. 
|4]and|5| [21 . The general trend is, both with and without H2O, that CO binds 
significantly stronger on top of Au2-4/MgO than oxygen does. 

As we showed in the previous section, the most favorable adsorption sites 
for H2O are in the proximity of the gold clusters, on top of the next-nearest Mg 
atom. The solid squares in Figs. [Ijandjsjshow the Eads [CO] and Eads[02] on top 
of Aui_4/MgO with H2O in the proximity. Some structures are energetically 
indifferent to the proximity of H2O, e.g. Aui^4C0 and AU2-4O2. In other cases 
CO or O2 bind stronger to the gold, e.g. AU3CO and AU1O2. Only in the case 
of the dimer CO binds weaker with H2O in the proximity. Figures |6] illustrates 
the ground state structures of CO/O2 adsorbed on Aui_4/MgO with II2O in 
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the proximity. The water molecule leads only in three of the systems to notable 
structural changes. First, the breaking of a Au-O^ bond in favor of a Au-H 
bond for AusCO/MgO, which is accompanied by charge redistribution within 
the cluster. Second, Aui^302/MgO bend down towards H2O and form a H-0 
bond. In Au302/MgO it also leads to a charge redistribution within the cluster, 

1. e. one of the Au atoms close to the surface loses its additional charge. 

Although the adsorption sites of H2O on top of the gold structures are ener- 
getically less favorable than on the surface. Fig. |3](a), we study the possibility 
for CO and O2 to co-adsorb on top of Aui_4H2 0/MgO (solid triangles in Figs. 
|4]and[5|. Carbon monoxide adsorbs on Aui^3H20/MgO as weakly as on the 
bare MgO surface. On the even-numbered clusters CO can bind stronger than 
on MgO (AU2) or essentially as strong as in the absence of H2O (AU4). The 
tetramer has two adsorption sites, the two gold atoms at the top. If CO and 
H2O co-adsorb on one of these sites each, they barely influence each other, as 
it can be seen, for instance, from their unaltered adsorption energies. 

The oxygen molecule can only weakly co-adsorb on Au2ll20/MgO, which 
does not lead to any structural changes. In the other three cases O2 forms, 
to different degrees, new molecular complexes with the pre-adsorbed H2O. Es- 
pecially on the odd-numbered clusters hydroperoxyl-hydroxyl complexes, i.e. 
02H--0H, are formed on top of the gold clusters. Fig. [7| by moving a hydrogen 
atom towards the oxygen molecule. Although the co-adsorption of O2 and H2O 
on the two adsorption sites of Au4 has a synergetic effect on their Eads, see Fig. 
|5j a hydroperoxyl-hydroxyl complex, as on Aui^s, is not formed. The hydrogen 
atom remains in the II2O molecule with a bond length stretched by 3% and also 
the intramolecular bond in O2 is only stretched by 6%. 

To complete this study of co- adsorption, we also address the question of II2O 
co-adsorbing on pre-adsorbed CO/O2 (open triangles in Figs. [4]and[5|. From 
the adsorption energies obtained for these scenarios one sees that H2O has a 
meta-stable local minimum on top of Au2,3CO/MgO and binds weakly on top 
of Au2,302/MgO, in the latter case without a significant additional activation of 

02. On AU1CO/O2 the water molecule can bind, although without a significant 
additional activation of the molecules. Co-adsorbed on the second adsorption 
site on top of Au4 we found the highest Eads[ Il20] with -0.64 eV. 

In conclusion, we predict the formation of the following ground state struc- 
tures under a mixed atmosphere: Aui02/MgO and Au2-4CO/MgO with a wa- 
ter molecule binding on the surface in the proximity of the cluster, cf. Figs. [3] to 
|5j Except for the dimer the adsorption energies of CO and O2 are enhanced by 
the proximity of a water molecule. Since the preferential gold/CO and gold/02 
structures are similar to those found for the water free case [21 , we conclude 
that also the blocking of Aui,2/MgO by O2 and CO, respectively, as well as the 
CO/O2 co-adsorption on Au3,4/MgO are essentially unaltered. Therefore also 
the catalytic activities of these systems are indifferent to the proximity of H2O. 
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5. Is a H2O mediated CO oxidation reaction possible on Aui_4/MgO? 

In Ref. [20 a H2O mediated reaction pathway towards CO oxidation has 
been suggested for Aug/MgO. The proposed catalytic cycle consists of the fol- 
lowing steps: 

Aun/MgO + H2O + 02^ Aun02H • • OH/MgO (1) 

Aun02H • • OH/MgO + CO ^ Aun(OH)2/MgO + CO2 (2) 
Aun(OH)2/MgO + CO ^ Aun/MgO + H2O + CO2. (3) 

It essentially depends on the formation of a highly activated complex, con- 
sisting of a hydroperoxyl and a hydroxy 1 group, i.e. O2II • • OH, on top of the 
gold cluster. Although we found that these complexes are stable on top of 
Aui^s/MgO, see Eads in Fig. [sj their formation probability is limited by the 
fact that H2O has only a meta-stable adsorption site on Aui^a/MgO, Fig. [3](a). 

Figure [7| shows the ground state structure of these two complexes and the 
calculated Bader charges. One hydrogen atom is moved from the pre-adsorbed 
water molecule towards the oxygen molecule, which redistributes its electron 
towards O2. The so formed hydroperoxyl group has a highly activated 0-0 
bond, i.e. the bond is stretched by 9% to 1.34 A. 

The remaining hydroxy 1 group receives an additional charge of 0.5 e~ from 
the gold clusters, i.e. a Bader charge of 7.5 e~ on OH instead of 7 e~. Fig. [t] 

Naturally the question arises, if these complexes open a reaction pathway 
towards CO oxidation as in the case of Aus/MgO? To answer this, we calculated 
the total energies of the system during the reaction steps of the catalytic cycle 
Eqs to (|3|, see Figs. [8]and|9l 

For Aui^s/MgO, where the O2H • - OH can be formed, we find it unlikely 
that a full catalytic cycle can be established. The formation of the first CO2 
molecule from a CO molecule and an oxygen atom from the hydroperoxyl group 
involves an activation barrier. We obtained activation energies for this reaction 
of 0.5 eV for Aui and 0.6 eV for Aua from a restricted molecular dynamics 
simulation. The kinetic energy of the CO molecule has to exceed these activation 
energies. Thus, a necessary condition for the oxidation reaction to happen is 
that the strength of the bonds B2 - B4 in the molecular complex, shown in the 
insets of Figs |8] and |9) exceeds the kinetic energy of CO. This might occur for 
Aui/MgO, since the remaining bonds B2 - B4 in the O2H • • OH complex have 
bond dissociation energies of at least 1.27 eV. However, the remaining complex 
(OH) 2 on Aui/MgO undergoes a barrier free dissociation into two hydroxyl 
groups. One of which is adsorbed on the adatom and the other one nearby on 
the MgO surface. This dissociation is possible since bond B3 is weakened to 
-0.6 eV once the oxygen atom is removed from the hydroperoxyl group. On 
MgO, on the other hand, the adsorption energy of an OH group is -1.19 eV. 
Both OH are inert towards a reaction with the second CO molecule, therefore 
interrupting the catalytic cycle. 

In the case of Aus/MgO the situation is different. First of all the adsorption 
energy of the O2H • • OH complex is more than 50 % smaller than that on the 
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adatom. At the same time the activation barrier for the formation of CO2 is 
higher, making this step even shghtly endothermic, cf. Fig. [9j This reaction 
becomes even more problematic, if one considers the bond dissociation energies 
of B2 - B4, cf. insets in Fig. [9j which are very close to the activation barrier. 
Hence it is quite likely that one of the bonds B2 - B4 breaks upon the impact of 
a CO molecule with a kinetic energy, necessary to overcome the reaction barrier. 
Still, if we assume for a moment that the bonds B2 - B4 will hold and the first 
CO2 molecule is formed, than the remaining steps in the proposed catalytic 
cycle are exothermic and barrier free, see Fig. [9j 

In conclusion, the dissociation of (OH) 2 in the case of Aui and the probable 
breaking of one of the bonds B2 - B4 in the molecular complex on Aua makes 
the proposed water mediated catalytic cycle unlikely to occur on the small gold 
cluster considered here. 

6. Summary and conclusions 

By means of density functional theory we studied the adsorption of H2O 
on small gold clusters, Aui_4, supported by a regular MgO terrace. We found 
that even these small gold clusters attract water molecules to their proximity. 
Two sets of stable adsorption sites exist: first, the energetically more favorable 
site on MgO in the proximity of the gold clusters, and, second, on top of the 
clusters. The latter showed a pronounced odd-even oscillation in the binding 
mechanism between H2O and gold, depending on the number of gold atoms in 
the cluster. 

We systematically studied the co-adsorption of H2O with CO and O2 on 
Aui_4/MgO and found that it is in general energetically more favorable to have 
CO/O2 adsorbed on top of the cluster with H2O in the proximity. Although we 
found the adsorption energies of CO and O2 to be enhanced by the proximity of 
H2O, the known small catalytic activities of Aui_4/MgO towards CO oxidation 
remain unaltered in this configuration. 

On Aui^a/MgO the oxygen and water molecules can form a stable, highly 
activated O2H • • OH complex. In contrast to the previously published findings 
for Aug/MgO, these complexes are unlikely to open a catalytic reaction pathway 
towards CO oxidation. The reaction barrier for the first CO oxidation is 0.6 
eV on the gold trimer, which is comparable to the bond dissociation energies in 
three of the five bonds in the 02H- • OH complex. Hence, the molecular complex 
is likely to decompose as a result of the CO impact. 

Although the bonds are stronger in the O2H • • OH complex on the adatom, 
the intermediate (OH) 2 complex, formed by oxidizing one CO molecule, disso- 
ciates barrier free into two hydroxyl groups, one on the adatom and one on the 
surface nearby. Since these OH-groups are found to be inert towards CO, the 
catalytic cycle is interrupted. 

We conclude that H2O can promote the co- adsorption of CO and O2 on 
Aui_4/MgO, but a stable, H2O mediated, catalytic cycle towards CO oxidation 
is unlikely to occur on these clusters when supported by a defect free MgO 
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terrace. In real-world samples surface defects, like steps and vacancies, might be 
responsible for catalytic activities of sub-nanometer sized gold clusters. Further 
investigations into this field are desirable. 
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Figure 1: (Color online) Ground state geometries of Aui_4/MgO(100) with H2O adsorbed on 
the surface in the proximity of the cluster. Note that H2O prefers to bind to the next-nearest 
surface Mg. Also shown are the Bader charges of the gold and water-oxygen atoms (in e~). 
The Bader charge of a neutral gold atom is 11 and that of the oxygen atom in H2O is 8. The 
different atom species are illustrated as yellow (Au), red (O), green (Mg), and blue (H) balls, 
respectively. 




Figure 2: (Color online) Ground state geometries of Aui_4/MgO(100) with H2O adsorbed on 
top of the clusters. Note the different binding mechanisms depending on the odd-even number 
of gold atoms in the cluster. Also shown are the Bader charges of the gold and water-oxygen 
atoms (in e~), cf. Fig. IT] The different atom species are illustrated as yellow (Au), red (O), 
green (Mg), and blue (H) balls, respectively. 
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Figure 3: (Color online) Top panel (a): Adsorption energies of H2O in the proximity (solid 
squares) and on top of Aui_4/MgO (open circles). The dashed horizontal line shows Eads of 
H2O on bare MgO. Lower panel (b): the separation energy Esep = Eo[Aun/MgO with H2O in 
proximity] - Eo[Aun/MgO] - Eo[H2 0/MgO] (solid triangles). Note that these points coincide 
with the solid squares in panel (a), i.e. E^ds [Aui_4] is not affected by the proximity of H2O. 
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Figure 4: (Color online) Adsorption energies of CO on bare MgO (dotted horizontal line) and 
on Aui_4/MgO (open circles). In the presence of a water molecule different co-adsorption 
scenarios are possible: H2O pre-adsorbed in the proximity of Aui_4/MgO (solid squares), 
H2O pre-adsorbed on top of Aun/MgO (solid triangles), and H2O adsorbed on top of 
Aui_4CO/MgO (open triangles). 
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Figure 5: (Color online) Adsorption energies of O2 on bare MgO (dotted horizontal line) and 
on Aui_4/MgO (open circles). In the presence of a water molecule different co-adsorption 
scenarios are possible: H2O pre-adsorbed in the proximity of Aui_4/MgO (solid squares), 
H2O pre-adsorbed on top of Aun/MgO (solid triangles), and H2O adsorbed on top of 
Aui_402/MgO (open triangles). 




Figure 6: (Color online) Ground state structures of CO/O2 and H2O co-adsorbed on 
Aui_4/MgO. In the three cases were new bonds are formed also the Bader charges are shown 
(in e-), cf. Figs, [l] and [2] In the gas phase carbon monoxide has Bader charges of 2.25 
(C) and 7.75 (O) and the oxygen molecule of 6 per O-atom. The different atom species are 
illustrated as yellow (Au), red (O), black (C), green (Mg), and blue (H) balls, respectively. 
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Figure 7: (Color online) On top of Aui^s/MgO highly activated complexes, consisting of a 
hydroperoxyl and a hydroxyl group, i.e. O2H • • OH, can be formed from an oxygen molecule 
and a pre-adsorbed water molecule, cf. Fig. [s] The Bader charges are shown as well (in e~). 
The different atom species are illustrated as yellow (Au), red (O), green (Mg), and blue (H) 
balls, respectively. 
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Figure 8: (Color online) Change in the total energy during the successive adsorption of H2O, 
O2, and CO on Aui/MgO. After the formation of a CO2 the remaining complex break apart 
into two OH groups, one on the adatom and one directly on the surface nearby. The upper- 
right inset shows the bond dissociation energies of Bl - B4 (solid triangles) compared to the 
activation barrier (dotted line) of the O2H • • OH complex, formed in step 3 and illustrated in 
the lower-left corner. 



14 




Figure 9: (Color online) Change in the total energy during the successive adsorption of H2O, 
O2, and CO on Aus/MgO to model the full catalytic cycle. The inset in the upper-right 
corner shows the bond dissociation energies of Bl - B4 (solid triangles) compared to the 
activation barrier (dotted line) of the O2H • • OH complex, formed in step 3 and illustrated in 
the lower-left corner. 
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